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After the completion of data taking at HERA-2, a large data set is at hand to study the photoproduction of charm 
and beauty quarks in ep collisions. New measurements of charm production based on D* meson tagging and beauty 
production based on muon and electron reconstruction test perturbative QCD calculations with improved accuracy. 
In general, QCD calculations at NLO describe the data well. The scale uncertainties are, however, large and dominate 
over the experimental uncertainties, calling for more precise calculations. 



1. INTRODUCTION 



O 

The study of charm and beauty production is a central topic of research at HERA. In ep collisions, charm and 
beauty quarks are predominantly produced by the photon gluon fusion process, which makes the production rates 
directly sensitive to the gluon density in the proton. The mass of these heavy quarks provides a hard scale so that 
calculations in perturbative QCD are expected to be reliable. 

However, this seemingly simple picture of photon gluon fusion with a hard scale given by the quark mass is 
considerably complicated by two issues: on one hand by the presence of additional hard scales such as the quark's 
transverse momentum, which means that heavy flavour production at HERA is inherently a multi-scale problem; 
on the other hand by the hadronic structure of the photon that leads to additional contributions such as flavour 
excitation, which pose additional challenges to theoretical calculations. 

While it is possible to take the stance that the theoretical issues are sufficiently well under control that one can infer 
• from the measured charm and beauty production rates the gluon density of the proton, a complementary approach 
C""^ , is to use the gluon density as determined from other processes, i.e. fits to inclusive structure function data, in 
particular from HERA, and investigate the extent to which the formation of heavy quarks is theoretically understood 
in a regime where multiple scales and QCD dynamics play an important role. 

In the following, new measurements of charm and beauty production in photoproduction at HERA will be pre- 
sented. Photoproduction is characterized by a near zero virtuality 1 Q 2 « GeV 2 of the exchanged photon, which 
essentially behaves like a real photon. 



X: 

2. THEORETICAL MODELS 



To correct experimental data for detector effects, Monte Carlo generators based on leading order (O (a s )) matrix 
elements and matched parton showers are used. The Pythia program |l] is based on collinear factorization and 
DGLAP evolution of parton densities. In contrast, Cascade [2^ is based on the fcx factorization ansatz, which uses a 
&t unintegrated gluon density that is evolved according to the CCFM evolution. 

For calculations at next to leading order, the FMNR program [3j, based on the calculation by Nason, Dawson and 
Ellis [4], is available. In this program, which does not include parton showers, the fragmentation to hadrons is usually 
modelled using an independent fragmentation approach; a new development is the FMNR <§5 PYTHIA program [5(, 
which augments FMNR with full Lund string fragmentation. 



lr The experimental definition is usually Q 2 < 1 GeV 2 and is driven by the detector acceptance for scattered electrons. 
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Figure 1: The HI measurement of D* photoproduction, compared to two QCD models at NLO: The FMNR program in the 
FFNS (shaded area) and a new calculation in the GM-VFNS (hatched) . 



3. CHARM PRODUCTION 

The HI collaboration has recently released a new measurement Q of charm photoproduction based on an integrated 
luminosity of 93 pb _1 of data taken at HERA-2. Charm events were tagged by the presence of a D* meson decaying 
in the golden decay channel D* + — > D°ir + . D° — > K~ir + (and charge conjugates). The data were taken using a 
new trigger 0] that allows a full mass reconstruction of the D* candidates at the third trigger level within 100 [is. 
The data covers a kinematic range defined by photon virtualities of Q 2 < 2 GeV 2 and photon-proton centre of mass 
energies of 100 < W 7P < 285 GeV. The D* meson is measured for transverse momenta pt(D*) > 1.8 GeV and 
pseudorapidities |r?(D*)| < 1.5. Differential distributions in p^ (D*), ij (D*) and W 7P as well as double differential 
cross sections in px (D*) and 77 (D*) have been measured and compared to various QCD models. In particular, two 
models at NLO QCD are available, namely FMNR and a recent calculation in the GM-VFNS Both models show 
similar behavour (Fig. [I]): The px (D*) spectrum is well described, with a slight deficiency at high px in the case of 
the GM-VFNS model, with sizeable uncertainties from scale variations, in particular at low pt- When the rj(D*) 
spectrum is considered, the theoretical uncertainties are several times larger than the experimental ones; nevertheless, 
both models show a tendency to fall below the data at large values of pseudorapidity. A similar behavour is found 
in the deep-inelastic regime [9j. 



4. BEAUTY PRODUCTION 

The measurement of beauty photoproduction at HERA is difficult due to the fact that beauty production corre- 
sponds only to about 0.03 % of the total photoproduction cross section. Two basic techniques have been used to 
enrich beauty production events: The measurement of track impact parameters or secondary vertices, using a silicon 
vertex detector, which enrich beauty production due to the relatively large lifetime of beauty hadrons, and measure- 
ments based on leptons (muons and/or electrons), which are based on the large semileptonic branching fraction of 
beauty hadrons. 

The HI and ZEUS collaborations have recently released results from several analyses based on lepton tags. Two 
analyses, one from HI [10( and one from ZEUS [111], are based on samples with at least one muon and two jets in 
the final state; in both cases, jet transverse momenta above p^ 1 ^ > 7(6) GeV are required. The resulting samples 
are enriched to about 30 % of beauty events by the requirement of two high-px jets in conjunction with a highly 
energetic muon. The actual fraction of beauty events is then determined from fits to the distribution of the impact 
parameter S of the muon track and the relative transverse momentum p™ 1 of the muon with respect to the nearest 
jet axis. In these measurements, the p™ 1 distribution allows the separation of beauty induced events from those of 
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light (uds) and charm quarks with good precision; the impact parameter distribution on the other hand allows in 
addition the separation of light and charm quark induced events and thus serves not only as a cross check of the p 1 ^ 1 
based measurement, but also helps to reduce the systematic uncertainty which would otherwise be incurred by the 
uncertainty of the charm quark contribution to the pJj? 1 spectrum. Both measurements cover a range of 0.2 < y < 0.8 
in inelasticity, momentum transfer Q 2 < 1 GeV 2 , muon transverse momentum p!^ > 2.5 GeV and jet pseudorapidity 
| 7 yjeti(2)| < 2.5; however, the ZEUS data cover a range —1.6 < rf- < 2.3 of muon pseudorapidities, while the HI 
analysis is restricted to —0.55 < 77^ < 1.1. The ZEUS analysis is based on 124 pb -1 of data taken in 2004 and 2005, 
while the HI result uses 171 pb -1 of data from the years 2006 and 2007. 

Both analyses show good agreement with perturbative QCD calculations performed with the FMNR program, after 
corrections for hadronisation effects have been applied to the FMNR predictions, which are available only at parton 
level. The excess of data over NLO predictions at low values of jet and muon transverse momenta that was reported 
in an earlier HI analysis of HERA-1 data and was in disagreement with ZEUS data [ijj], is not confirmed by the 
new HI analysis. 



ZEUS has also published a new measurement based on events with semileptonic decays to electrons 14|. In this 



analysis, events with two jets of E^ 1 ^ 2 ' > 7(6) GeV and at least one electron are selected from 120pb _i of data 
taken at HERA-1. The beauty contribution was determined from a fit to a likelihood variable based on quantities 
that are sensitive to the quality of electron identification, such as the specific ionisation AE /Ax, as well as quantities 
that allow the separation between electrons from light quark or charm induced processes from bottom production. 
These quantities are again , and Acf>, the azimuthal angle difference between the electron flight direction and the 
direction of the total missing transverse momentum. Again, reasonable agreement with the NLO QCD calculation 
FMNR is observed, for the total cross section as well as differential distributions. 



An upcoming publication by ZEUS [15| is based on 
114 pb -1 of data taken at HERA-1 and utilizes events 
with two muons in the final state. Contrary to other 
work discussed here, no jet requirement is imposed 
in the selection, which is possible due to the high 
beauty fraction of a dilepton sample, and the anal- 
ysis is performed down to muon transverse momenta 
of pt^ > 1.5 GeV, or even pt^ > 0.75 GeV for high- 
quality muon candidates. The analysis is thus sensi- 
tive down to to the threshold of beauty production, 
where the transverse momentum of the beauty quarks 
is close to zero. Therefore a measurement of the to- 
tal beauty production cross section is possible with 
relatively small extrapolation errors. The total cross 
section for the process ep — ► bbX at y/s = 318 GeV 
has been determined to be a to t(ep — ► bbX) = 13.9 ± 
1.5 (stat.)!^ (syst.) nb, to be compared to the NLO 
QCD prediction of <r^°(ep -> bbX) = 7.5 ±±£f nb. 
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Figure 2: A compilation of beauty photoproduction cross sec- 
tion measurements as a function of the mean transverse mo- 
mentum of the beauty quarks, compared to QCD calculations 
at NLO. 



Within the large uncertainties, in particular of the 
NLO calculation, the NLO prediction is consistent 
with the data. 

Fig. [2] summarizes all currently available beauty photoproduction measurements at HERA, in comparison to the 
QCD prediction at NLO. It can be seen that NLO QCD describes the data well in shape and normalization, albeit 
with rather large uncertainties, which are dominated by the choice of scale. Note that for this plot the central value 

i/2V^T~ 



of the QCD prediction is derived for a choice of (j, t j = /iq 



Pt 



scales [l6|, while the individual measurements sometimes use twice this value, i.e. fio — J m? + p\ for the central 



value. 
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5. CONCLUSIONS 

The study of the photoproduction of charm and beauty quarks at HERA remains a challenging testing ground for 
perturbative QCD. With the increased statistical precision provided by the HERA-2 data set, and improved detector 
understanding that results in more precise measurements from HERA-1, measurements of charm and beauty produc- 
tion are now available with experimental uncertainties that are significantly smaller than the theoretical uncertainties 
of QCD calculations, even at next to leading order. Using gluon densities based on fits to inclusive structure function 
measurements, these calculations in general describe the data reasonably well in shape and normalization, which 
is a great success of perturbative QCD. However, differential distributions such as the n distribution of D* mesons 
indicate deficiencies in the calculations which may need further work to resolve. 
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